General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



NASA TECHNICAL 
MEMORANDUM 


NASA 7 * 



Report No. 53909 


PRODUCT MEAN VALUES AND CONVECTION SPEED 


By W. H. Heybey 
Aero-Astrodynamics Laboratory 


September 12, 1969 





NASA 

George C. Marshall Space Flight Center 
Marshall Space Flight Center ; Alabama 

(THRU) 

zfei 

(CATEGORY) 


*70-375,65 

S (^CESSibN^^B ERT^ 


o 

<o 

Q£ 

O 


MSFC - Form 3190 (Stpt.mbtr 1968) 


(PAGES) 

.. TMzIW-, 

6 (SaSACRORTwS ORAD NUMBER) 


NASA - GEORGE C. MARSHALL SPACE FLIGHT CENTER 

changed to TM X-53809 on September 12, 19£9 


IN-AERO- 67 -8 


November 20, 1967 


PRODUCT MEAN VALUES AND CONVECTION SPEED 

By 

W. H. Heybey 


This note has been prepared for internal use within the Aero- 
As trodynamics Laboratory and is transmitted to other agencies only for 
expediency in special cases. Transfer to parties other than those 
approved by Aero-As trodynamics Laboratory Directory is not authorized. 


TECHNICAL AND SCIENTIFIC STAFF 
AFRO- AS TRODYNAMIC S LABORATORY 
RESEARCH AND DEVELOPMENT OPERATIONS 


a»mi 


IN-AERO-67-8 


PRODUCT MEAN VALUES AND CONVECTION SPEED 


By 


W. H. Heybey 

George C. Marshall Space Flight Center 
Huntsville, Alabama 


ABSTRACT 


In crossed-beam methodology the relationship of the beam intensity 
covariances and local extinction coefficient covariances comes up as an 
important problem. The approach leading to it has been carefully 
restated in a systematic review. The answer suggested places more 
weight on experimental determination of correlation volumes and flow 
direction than on mathematical conditions imposed a priori on the 
turbulent state. Application is made to the calculation of convection 
velocity. 
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PRODUCT MEAN VALUES AND CONVECTION SPEED 

SUMMARY 


The physical conditions under which basic mathematics currently in 
use are valid have been carefully specified. If they are realized, 
crossed-beam examination of a thin turbulent layer promises to yield 
reliable results. In thick layers the turbulent correlation must remain 
restricted to small volumes if one wants to calculate two-point values 
from two-beam values. The experimentation normally necessary to explore 
the flow's correlation structure can be forgone in computing bulk veloc- 
ity, if one agrees to define it by the time of the maximum two-beam value 
and the beam distance in flow direction. The latter must be ascertained 
as (approximately) constant simultaneously with any bulk speed measure- 
ment . 


INTRODUCTION 


The relationship of the two-beam and two-point product mean values 
is considered a basic issue in crossed-beam experimentation. Since it 
has been controversial, and still is, to some extent, an attempt is made . 
to examine the case in a systematic way, starting out with the very 
foundations of the crossed-beam method in order to ascertain the physical 
significance of the mathematical formulations that follow from them. 
Determination of the convection' speed is always kept in view, as it is a 
main application. 


« 

ORIENTATION 


Two non-parallel light beams, a and b, of infinitesimal width are 
sent through turbulent flow at some distance from each other. Let us 
mark on them the points A and B as the origins of running coordinates, 
a and {3, defining the points on the a- and b-beams. The flow state 
along the beams will have to be described by functions that depend on 
either d or (3 and on time. 
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It will be assumed from the outset that the convection velocity has 
a constant direction given by the unit vector v . If it is not parallel 
to the pair of parallel planes that is established by the lines a and b 
(as we will postulate), there exists exactly one pair of points, P A and 
and P s , on the beams such that P A Pg ;j v. Since these points are in line 
with the flow, the physical correlation near them will be stronger than 
near any other pair. We select them as the origins A and B, both of 
which then have a definite location. It follows that their connecting 
line will not in general be normal to a or to b. 

We wish to inquire into what can be learned through the modifica- 
tions of the beam intensities alone, and are not therefore concerned 
with the state of the flow outside the beams, To view the a- and b-lines 
as parallel to two axes of some rectilinear three-dimensional system is 
not necessary as a consequence, nor has it been found advisable to do so. 
Especially, the line AB = s is not considered as parallel to or coinciding 
with a third axis. While the distance s is constant for every experi- 
mental run, a and £ are true variables, as the light is acted upon at 
every point P A , and at every point P B . 

A beam, operating for a time span of T 0 seconds, can be imagined as 
an infinitely dense sequence of separate instantaneous "f lashes." A 
flash "at the time t," needing practically no time to travel from the 
source to a detector, will nevertheless be weakened through various 
absorption and scattering mechanisms that it encounters on its way, e.g., 
along the a-trace. Its intensity, I a , recorded M at time t," varies with 
t, since in a turbulent flow those mechanisms work in a slightly differ- 
ent manner on succeeding flashes. The detector record will exhibit the 
function I a (t) as a random curve over time. The temporal mean, taken at 
time t*, will, with sufficiently persistent turbulent characteristics, 
become practically independent of t* after t* has grown sufficiently 
large, say, from t* = T* on. Mathematically, this type of turbulence 
is described by what may be termed the "persistency" condition: 

T 

I a = - jT l a (t) dt pb const. for T ^ T* . (1) 

o 


Failure of a record to comply \-7ith this condition may be caused either 
by too short an operating time* (which deficiency might be amended) or 
by macroscopic time trends in the turbulent state for which there is no 
remedy. 


It is clear that T cannot be taken larger than T 0 . 
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It may be mentioned in passing that, if a whole series of a-records 
of lengths T 0 were taken in rapid succession, there is no guaranty that 
the ensemble means should be the same at all times (i.e., stationary), 
merely because, in every sample, condition (1) is satisfied. 

The absorptive weakening of light intensity is described in optics 
by the exponential law 


dl - -K I da. 


( 2 ) 


In a homogeneous gas (as considered in optics) the absorption coefficient 
K is a constant for a given wavelength and refractive index (density). 

The latter varies in steady (but not uniform) motion, so that the gas 
is "homogeneous" in its smallest parts only; in turbulent motion, it is 
so only instantaneously. While maintaining the law (2) and settling for 
a definite wavelength*, we must here regard K as a function both of a 
and the time t at which the flash occurs. The latter’s intensity, I. 
along the a-beam likewise is not only a function of flash time, but also 
of location, I = I (a, t) , If a^. and a 2 are the coordinates of tie sender 
and detector, respectively, I(a l3 t) = I s (t) , I(a 2 >t) = I a (t)«. The first 
formula allows for possible time fluctuations of the radiated intensity. 
However, for the purpose of the present theoretical investigation, ideal 
operation will be presupposed: there are no power fluctuations at the 

source, no instrumentation or other noise effects. Thus, I g (t) = I s = 
constant. 

The law (2) will be assumed to hold as well for the iocal intensity 
loss caused by any other dissipation process, as by molecular** or 
particle scattering; K will therefore be given the more general name of 
"extinction coefficient." 


One can use white light as well. This Is not done in optics, because 
there the interest precisely centers about the questions which fre- 
quencies are absorbed and to what degree. 


** Lord Rayleigh has found that here 


K 


.( s . a - 1 ): 

NA 4 


(n - refractive index, N '= number density, A - wavelength). 
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On integrating the differential (2) 


I (t) 

d 


a- 


exp 


K(a,t) da T . 


_ i.' | I \'"“J ~t — VOV I 

I S l J j 

a x 


(3) 


The average of I a over the time span of the integral (1) follows as 


a 2 T 
nr. 


0 fe 


\ = ^ eKp 


L T 


K(a, t) dt 


Q$i o 


da = 


= I exp 


K(a) dar , (4) 




where K(a) is the time average of K at station a. I could be taken here 
as equal to the operating time T 0 ; but it is not advisable to do so. 

With the b-beam, a time delay t will be contemplated, while the integra- 
tion still has to extend from t = 0 to t = T. 

If the difference 


k(a, t) = K(a, t) - K(a) 


(5) 


is introduced, the ratio (3) becomes 

1 (t) r I r F 2 . 1 

— — — = exp -j- / k(a, t) daj- • exp -j- J K(a) daj- , 
s a! a x 

and, with the use of relation (4), shifts into the ratio 


I (t) 


exp 


a 2 (Xq _2 

■ jT k(a, t) daj- =1 - jT k da + ^7* J k da 
Qh . <^1 


-f- 

1 • • f • 


( 6 ) 


Suppose now that the turbulent motion is such that the assumption 




k(a, t) da 


« 1 




(7) 


can be made for every time t. Then it follows from the ratio (6) that 


a 2 

i(t) a I - i (t) = x / k(a, t) da. (8) 

a a a J 


It is seen that the linearizing condition (7) implies that the recorded 
intensity I a (t) should never stray far from its mean value. This does 
not necessarily exclude rather large values of k(a,t), that is, marked 
deviations of the local value K(a, t) from its temporal mean K(a)> since 
for every flash they might effectively cancel out when integrated over 
O'. But, whether the a-beam cuts through vigorously fluctuating or rela- 
lively calm flow, if large-scale temporal trends in the extinction pro- 
cesses are present at many or all points P^_, the detector record will 
reveal that I&(t) is not close to I a everywhere* Linearization is then 
dubious, if not prohibited. The persistency condition (1) which is likely 
to be violated in such records with long-term up-and-down or monotonic 
trends, can be expected, on the other hand, to be satisfied if but small 
deviations occur from the intensity mean. We may say that if the persist- 
ency condition is observed as realized, the linearizing assumption (7) may 
be taken for granted. Even if the record fails to comply with this condi- 
tion, linearization still may involve no risks when long-term temporal 
trends are not apparent. (Those trends are sometimes classified as non- 
‘s tat ionary -- another use of an overburdened word.) 

As an immediate consequence of definition (5), 


k(a, t) dt = 0, 


(9) 


Both the requirements (7) and (9) should be heeded when setting up 
approximating expressions for the function k(a, t). It is already 
recognized that T should as a rule be smaller than T 0 , The limits of 
the integral (7) can be narrowed down in circumstances. It is not neces- 
sary that the turbulent realm should extend all the way from sender to 
detector, provided that outside it the extinction coefficient is con- 
stant with respect to time. Let it be equal to Kx(Oi) in % < a a 
to K‘-'( a ) i n a < cq> Then 


V 0 


a\ 

/*v *** 


Go 




exp -j- / Ki(a) da - / K(a,t) do; - / K L >(a) daj 


«i 


r 

«1 


a’ 


and 


I 


a 


exp 


{-/ 

CC X 








KiCa) da - / K(a) dt - / K 2 (a) da, 




a\ 


that 


a 2 


exp 


L (t) 


’ P “* 

-J K(a,t) da ^ 


a 


P'2 


^2 -v 

exp K(a) da] 


- l - J k(a, t) da, 
a!, 




in accordance with expression (6). Definition (5) and condition (7) 
now refer to the turbulent interval < a^, a 2 > alone. In the future, 
the integration limits a 1? a 2 will be regarded as the bounds of the 
turbulent segment on the a-beam. The values K x and K s in the laminar 
parts are of no consequence, provided that they truly depend on O alone: 
The outside 'flow must be steady; number and size of scattering particles 
are not allowed to vary with time along the laminar segments. 
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Analogous expressions evolve with the beam b under the conditions 
stated. In writing them down, primes will be used if necessary for 
clarity. Thus 


T 

| Ji b (t + t) dt, 
o 


but 


F 3 

i ' (t + r) - I b ( t) - I b (t + t) = i b (T) J k 1 O, t t) dp, etc. 

Pi 


The crossing of the beams is performed, not indeed physically in 
general, but mathematically by forming the temporal covariance of 
I a (t) and X^(t + t); the value of i at the station t of the a-record 
is coupled with the value of i 1 at the station t + t of the downstream 
h-record. The delay time t is needed for later applications and has the 
character of a parameter on which the covariance (or "two-beam product 
mean value") depends; 


T 

n 


r (t) - ~ j i(t) i * (t + t) dt 


o 


s. ~ / \ T 0^2 

% p r p 

7jT~ — J J k(QJ, t) da 

o ax 


V (p, t + T) dp 


dt = 


Pi 


= i a ^(t) R(t) . 


( 10 ) 
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When the two records are laid above each other with coinciding time 
scales and time scale zeros, it is understood in this formulation that 
the b-record is used from t~ t to t ~ T + ;, with a number of parametric 
.-values. These of course have an upper bound t u * T 0 - T, while the 
lower bound is r = 0.* 


The function R(t) introduced above may be written as 


R(t) 



«1 01 


k(a,t) k' O, 


t + t) dt 


d£3da . 


( 11 ) 


The rearrangement of the integration sequence is permitted, since both 
the integrands and the integration limits are finite. One recalls that 
condition (7) and the analogous condition for k' ({3, t + t) require that 
the integrand (10) is numerically very small. 

The functions k(a, t) and It' (p, t + t) describe, relative to a tem- 
poral mean, the time fluctuations of the extinction coefficients at the 
two points P A (a) and Pg(f3), Accordingly, one may define temporal two- 
point product mean values 


0; T) = ~ J k(a,t) k' (p, t + t) dt, 


( 12 ) 


so that the function R(t) may be set into the form 

R k (a,(3;T) dpdoj. (13) 

Pi 

R| c is the covariance of K(t) and K* (t + t) calculated at the points 
P A (a) and P B (p), If the choice of these points is inconsequential, 


•» T - . 

“Negative values of r must be used when the a-record rather than the 
b-record is being "delayed. H 



if Rj c docs not actually depend on a and t ; ., so that 


i.e 


> 


R k (a, (3 ; t) - Rj (0,0; t) - const., 
we have the simple and desirable relationship 


R(t) = (a 2 - ai)O a ** $ 2 ) 1^(0, 0;t), 


(14) 


which shows the two- point mean values that would be found at A and B 
(or at any other two points P^ and P^) as proportional to the observable 
values R(t). That turbulence should ever be so organized in toto appears 
as highly improbable. It is, however, one prime aim of the crossed-beam 
method to gain local information in spite of the use of integrating 
light flashes. The question then arises: In what circumstances is it 

possible to relate the two-beam values R(t) to the unknown two-point 
values Rj c ^0,0;t) by means of r-independent factors of proportionality? 


SHEET TURBULENCE 


By definition, both the points A and B (and their connecting line 
as well) are inside the turbulence region. Let us assume that by 
chance they occupy the centers of the respective turbulent stretches, 
whose boundaries then may be given as 


» -Q*, a 2 = 05*; Pi = -p*. 


Pa “ P*. 


(15) 


This simplifies the results without blotting out essential features. 
The distribution at flash time t of the k- values can be approximated 
by the trigonometric polynomial 


M 


^ C a m(t) cos moi 4* G^(t) sin ma] . 


m=0 


k(a, t) = 


(16) 


If M is a sufficiently large number, one may presume that the polynomial 
will yield a satisfactory approximation to the distribution at all times 
t. The values of the fluctuations a m (t) and C%,(t) theoretically can be 
determined for a given turbulence field by running the gamut of all t's 
and finding the best approximation for every single flash along the 
a- trace . 

We postulate that k(a, t) satisfies the conditions (7) and (9). 

In the same vein we define 

k' O, t + t) = ) [b (t + T ) cos n P + f3 (t + t) sin n(3] . (17) 

n=0 


From these definitions we obtain the expression 






b (t + t) dt. 
n ' 


(18) 


On the other hand, with the limits (15) the function R(t) becomes 


R(T) 


4 

T 



sin mcT 
m 


o m 



b (t + t) 
n 


sin nB ,f 
n 


(19) 


There is obviously no constant factor of proportionality that would 
connect R k (0,0;T) to R(t) for all values of t except when both Of and 
P' r are very small. Then 


R(t) = tetf R k (0,0;ir). (20) 

The factor 4cTf3' f is the product of the turbulent segments along the a- 
and b-beams. One can show that, if the origins A and B are not their 
center points, the result (20) goes into the relation (14) when again 
the segments are very shprt. It does not matter then, where A and B 
are located, R k (0,0;ir) is always equal to R(t) divided by the product 
of the turbulent segments. (With approximations for k and k' other than 
by trigonometric polynomials the divisor may turn out merely propor- 
tional to this product.) 
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For defining dimensionless (normalized) expressions let us introduce 
a third beam, c, into the experimentation which physically crosses the 
a -beam at the point A (doubling as origin of a linear y-coord inate) , and 
which often (though not necessarily so) will be chosen as parallel to 
the beam b. If the a- and c-records are correlated without time delay 
we may introduce a fluctuation 

P i 

k" (7>t) = ^ [c p (t) cos py + 7 p(t) sin py] . 

p=0 


Note that even with c jj b and P - N the coefficients c^(t), y^(t) cannot 
in general be supposed to be equal to b^(t + 0), p^(t + 0). The dis- 
tance s being zero with the a- and c-beams, we will use a corresponding 
subscript in writing down the analogues to the formulas (18), (19), 
and (20): 


\,o^ a = °> 7 = 0 ; t = 0 ) = “ 


a (t) • ) c (t) dt 

, m w p w 


o m 


R q (t = 0) 


4 

T 


a (t) 
m 


sm m g' 
m 




w 


c (t) S2&2Z- dt 
P P 


o m 


In the event that c<? f and y* are very small 


R o (0) = 4aV f \ }O (0,0;0). 


The normalized values of R^(0,0 ;t) and R(t) will be defined as 


\(0,0;t) 

q = and 


Q 


UiL. = 

R o (°) 


i b (T) 


A/ 

R(t) 


( 21 ) 


( 22 ) 


so that, if a 5 % 3 " , 7 * are very small 



c 



( 23 ) 


In the light of the foregoing the method promises adequate results 
when diligently applied to thin turbulent sheets, as for example to 
thin free jet boundary layers anchored at the mouth of wind tunnel 
nozzles. The outside flow is steady, while the turbulent state can be 
expected to possess the overall persistency required to satisfy condi- 
tion (1) and thus to permit the linearized approach. Furthermore, 
close to the nozzle exit the direction v can generally be taken as 
parallel to the known (and constant) direction of the undisturbed flow. 
The form of the factors of proportionality appearing in relations (20) 
and (23) depends on the functions chosen to approximate k and k 1 ; how- 
ever, the mere fact that proportionality exists suffices in all cases 
where scale effects are irrelevant, e.g., in determining convection 
speed, which calls for finding that particular spot, t = t*, at which 
R| c (0,0,t), and therefore R(t), become maximal. Strictly speaking one 
must insure that I]-, (t) does not vary with t to any significant degree, 
since by expression ( 10 ) and ( 20 ) 

R (0,0;t) = — 7-7 — R(t). 

* toV i a i b (T) 


The convection velocity 




v 


is usually found somewhat smaller in magnitude than the undisturbed 
velocity parallel to it. 


VOLUME TURBULENCE 


The proportionality (14), if proven valid, relates the two-point 

value for A,B to the two-beam value along a,b. The important role of 

the direction v is to define those two especial points on the a- and 
*<>— » • • • • 

b -beams for which the proportionality has a chance to exist. The 
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strong correlation i c twee a tiu- origins will gradually subside as an* 
ii } | incrcuan . Outr.ido certain "correlation segments" surrounding A and B 
the contributions of R k to tin* value R(i) as given by expression (13) 

'viT < l e negligible. licit t g'.enis can be expected to be relatively 
narrow if a and b are far from parallel. Still, with unruly turbulence 
there will be sizable contributions to the integrand (13) by points out- 
side u « 0 and 3 = 0 with values R k (a, (3; t) R^(0,0;t) so that the pro- 
portionality (14) is precluded. 

To gain access to the problem posed therewith let us pass through A 
a number of a -beams , a( n ), in different directions and equally many 
b-beams, b ( n ) , through B, in directions more or less vertical to the 
plane containing AB and the corresponding beam a( n ). The two-point 
product mean values R£ n ' (a, ft; t) will usually be different for the 
several n's except at the points A and B which are common to every pair 
of beams. Setting 


; (0,0;t) = R k (0,0; T ) J 


we have the Taylor expansions 


R^ n) (a, p; t) 


\(0,0; T ) + « —§5- + P c^q + • • • 

3^=0 f^o 


(24) 


where the first term is the same for all pairs. 

/ \ 

Also, the covariance functions R' - } (t) differ in general from each 
other, as do the functions 

~(n) , s 1 „ (n) . . 

R (T) " =(n) -(n) R (T) - 
a Ti 


(It is again assumed that the means iy ' (t) do not noticeably vary with 

T.) 
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If, however, the experiments indicate that, for all n, the values 
of R( n > (t) are, al though t perhaps not equal to, so at least proportional 
to those of a function R“(:): 


R (n) (T) = K (n) r*(t) (25) 

it will also hold that 

R (n) (t) = A (n) R*(t). (26) 

It follows from expressions (13) and (24) that 


R (n) (T) 



2 _ 


(n) 




bp 


(n) 


dpda 


In order to establish relation (26) the right side must not depend on 

— (nor on the higher derivatives), meaning that ql and (3 
must remain very small. In other words, the integral (13) must be 
zero outside two very small correlation segments (d^ - a4 n ^), 

(p^ " Pi^), inside which is practically constant with respect to 
d and (3. The proportionality (26) then assumes the form 

R (n) (t) - (c 4 n) - ai n) )(Ps n) - pi n) )Rj t (0,0;T). 

so that the relation (25), observed as true experimentally, may now be 
written as 

R (n) (T) = n (n) R k (0,0,T). (27) 
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This is the result desired. Since it holds for all conceivable (a,b)~ 
combinations (for all n), there exist two minute correlation volumes 
surrounding A and B such that, at any given t, the Revalue is practically 
the same no matter which pair of inside points on two such beams is being 
correlated. The turbulence is then sometimes said to be isotropic (in 
the neighborhood of the pair A,B)*. 

The beam directions have a bearing on the factor p/ n ^only and can 
be chosen at will in velocity calculations since, by hypothesis, the 
maxima of all the functions R( n )('T) reside at the same value, t= t*. 
Indeed, if one is interested in this value only, the relations (25) may 
be required to hold merely in a narrow t*- environment. The proportion- 
alities (27) would then be true for values close to r* only, but that 
is all that is necessary for the purpose intended. 

One may argue that, with u constant bulk velocity direction, it 
might suffice to establish isotropy normal to it only, and accordingly 
may restrict the beams a( n ) and b ( n ) to those perpendicular to v (and 
to each other). In the event relations (25) are found to be true, the 
correlation "discs" surrounding A and B will have small diameters in all 
directions . 

Whether or not the experimentation described in the foregoing can 
be carried out in practice, will have to be decided by those who would 
be in charge of it. In any case some experimental means must be con- 
trived to make sure that in a turbulent flow where correlation volumes 
(or discs) might be sizable, they are in fact not so. 

Further experimentation is called for, should the direction v not 
be known beforehand. 

For example, one could operate with two parallel beams (not too 
close to each other) and determine the value R(0), which will be small, 
unless the plane of the beams is parallel to the bulk velocity when 
R(0) will have increased to a maximum. The plane's proper attitude 
can be found by rotating the a-beam about the axis b. A second pair, 
with a different direction of b, will secure a second plane, also 
parallel to flow direction; the line of intersection yields the direc- 
tion ± v. Repeated experimentation will show whether or not v meets 
with the requirement of constancy in the turbulent zone investigated. 


’'The usual definition of isotropy refers to one single point and to the 
time-averaged squares of the velocity fluctuation components in any one 
d irection. 
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Any plane known to be parallel to v, but not to the first plane, 
would alleviate the work load, since it can serve a'j the second plane. 
In the study of atmospheric turbulence it is often justified to take 
such a plane as horizontal. 


CONCLUDING REMARKS 


Physically, the requirement of small correlation volumes would 
mean that, if A and B are considered as the centers of eddies, their 
diameters must be very small. This is hardly surprising; beam correla- 
tion can be reduced to point correlation, when the correlating agents 
are point-like themselves. Moreover, the concept of convection speed 
has a clear meaning only then. It grows ill-defined with eddies of 
finite size, since every fluid particle residing in them pursues its 
own course, sharing a common or bulk velocity component with every 
other particle. One might be willing to define the latter with the 
use of that particular value, t = r*, that maximizes R(t), a function 
no longer proportional to R| c (0,0;t), but also no longer required to 
be so. An eddy particle leaving point A will not as a rule pass through 
B, nor even cross the line b. The geometric bulk of the eddy, however, 
will abide in the neighborhood of B after a travel time of about t* 
seconds. In these circumstances experimentation to prove isotropy is 
not necessary, while the direction vector v must be determined and 
found constant, so that the length AB ~ s is available for computing 
the bulk speed. 


